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I. Introduction 



Current interest in string theory in less than or equal to one dimension has arisen 
from the discovery of the existence of the continuum limit called double scaling limit 
]]J which sums string perturbation theory. At the same time, it has become clear 
that the universal equations which govern the correlation functions near the critical 
points are given by hierarchical integrable differential equations ^ f|. 

By some time last year, however, more than several people had come to realize 
that many of the results obtained in the limit are already visible when N- the size 
of the matrix- is kept finite. This provides an opportunity to study the system in its 
original definition and as a solvable model. The progress from that direction is what 
we will discuss below. We will summarize some aspects of zero dimensional matrix 
models from the point of view based directly on results obtained at finite N. 

Such analyses appear to be imperative in view of the recent status of nonpertur- 
bative 2d gravity: no limiting procedure has been found, up to now, which maintains 
both the reality of the partition function and the original combinatorial correspon- 
dence of a matrix model with a triangulated two dimensional surface ||, |5|. No sat- 
isfactory definition of continuum 2d gravity ( nonperturbative) has, therefore, been 
given. 

In the one matrix model, we will cover the following items: Virasoro constraints 
at finite N fully characterize the space of correlators: the equivalence of the two 
approaches, i.e. Dyson-Schwinger approach and the one based on orthogonal poly- 
nomials; how to take the double scaling limit of the Virasoro constraints at finite N 
near the (2, 2k — 1) critical points || to obtain the Virasoro constraints of a twisted 
boson [|7|, ||; correspondence between the one matrix model and the classical Toda 
lattice equation |1: the derivation of Kazakov's loop equation |IIJ directly from the 
Virasoro constraints (See, for instance, |Tl[| ). See also |L2] for these items. 

In the case of the two matrix model, the situation has been clarified since the 
original developments and proposals [0, |T5|| . By now, we agree that the two matrix 
models capture the essential feature of the 2d gravity coupled to a general (p, p') 
conformal matter |L4|, [15f1 : it contains all the critical points indexed by a set of 
coprime integers (p,p')- We will discuss the new constraints of type derived in 
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A. One matrix model 

Recall the partition function of the one matrix model: 
4 } KM}] = f /Me-^(^) 



A' 
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V(A) = £>A< . 

The Dyson-Schwinger equation is a set of consistency conditions on the space 
of correlators which is derived under general variations. It can be regarded as an 
equation for the space of string theories of a particular class ( in this case, the ones 
described by the one matrix model). An efficient way which leads to the Virasoro con- 

N d N 1 d 

straints is the following: insert — ^ K + ~JV~ or E 7 TTT anc ^ ex P ress the resu lts 

i= i d\i i=1 £ — Aj aXi 

in terms of the couplings and the derivatives of the couplings in two different ways, 
using partial integrations. We obtain, after some calculation on the Vandermonde 
determinant, 

e n z$ = , n>-l , 

oo Q n pp. 

L = ^ 9 ^ + ^d^k- ■ (2) 

These are the Virasoro constraints at finite N. In this derivation, it is evident 
that they arise from the reparametrization of the eigenvalue coordinates. 

As a natural extension of the above derivation and also as a warm-up to derive 
constraints of the two matrix model, consider the higher order differential operators 

13 

or Z^UY • (3) 



Let us introduce the following notations 
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In ref. ||16||, we succeeded in writing the constraints in a closed form: 
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The existence of such constraints is peculiar in view of the fact that the Virasoro 
constraints in the double scaling limit fully characterize a space of correlators of 
topological gravity |7], ||. We have checked |l6j in the lowest few cases that these 
higher constraints are in fact reducible to the Virasoro constraints: 



■u/ 2 ) (£) Virasoro 

«,<8>(0 =(jV 2 ) (0)^ + ^(0 , 
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In the approach based on the Dyson- Schwinger equation, we study a relationship 
of all possible singlet operators and their correlation functions. An alternative ap- 
proach is the one based on orthogonal polynomials, where one studies all the matrix 
( both diagonal and off diagonal) elements of the eigenvalue operator A as well as its 
conjugate 4^. 

Let us enumerate basic properties of the orthogonal plynomials P n (A) = A™ + 
c n -iA n 1 + ■ • • : 

XP n (X) = P n+1 + e^P n (X) + R n P n ^(X) , 

S n ,m = d\ e ( 9 ' — -=. -=- =<n\m> , 

h n +i — Rn+\h n , 
2V-1 

Z N = N\ J] h . (8) 
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The basic equations in this approach are the matrix elements of the Heisenberg al- 
gebra and a response of the system to a change of parameters, namely, a parametric 
derivative of the normalization constants h~: 
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(9) 
(10) 



where = (ft) are a set of parameters which are functions of the original (bare) 
couplings gi. The first equation (§) is regarded as a discrete version of the string 
equation. The second one (|I0| ) may be called a discrete flow equation. The basic 
equations derived in the continuum theory are already visible when N is kept finite. 

In order to demonstrate the equivalence of the two approaches, let us rederive the 
Virasoro constraints from the equations (pl),(|T0|) above. Begin with 



(n + 1)A T 



d_ 



A 



n+l 



N-l 



Take a trace of this relation and multiply by Z^. Namely, < j 
Using eq. (p70|) (tt = gi), we find 
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j > z N . 

(12) 



For the second term of the right hand side, use a formula which holds for any one-body 
operator 



j:<j\o(\,±)\j>=<<j:o(\ t ±) », 



i=l 



d\i 



(13) 



where << • • • >> ave denotes an averaging with respect to the partition function 
(eq. ([I])). ( The derivative is in between the two determinants and does not act on 
the potential.) Eq. (pID readily follows from 



A(A 1? • • • Aiv) = det (Pi-i(Aj)) 



(14) 



The calculation has reduced to the one mentioned in eq. @. We obtain the Virasoro 
constraints at finite N again. 
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The properties of the orthogonal polynomials have a curious connection || to the 
classical Toda lattice equation. Let 

4> n = lnh n . (15) 

From eqs. (|9|) , fllOD , we find a hierarchy of classical differential equations whose first 
member is the classical Toda lattice equation: 



( _JLV rfi = g^n + l-^" _ e <Pn + l-<Pn Qg\ 



We leave the details to the references 0. It is our hope that explicit results for the 
physical quantities come out from this correspondence at finite N. 

Note that we are still dealing with quantum theory. The feature that quantum 
correlators obey classical integrable differential equations is shared by other exactly 
solvable models of quantum field theory [IJJ . 



It is instructive that the Kazakov's original loop equation [10 can be directly 
derived from the Virasoro constraints without touching the properties of the measure 
factor i.e. the Vandermonde determinant. This fact has provided much inspiration 



to a recent work [IT] on the construction of the supersymmetric loop equation by 
Alvarez- Gaume, Manes and the author. Let a loop with length I be represented by 

JV 

*>(A = ]■■■ \T.f Xl ) e 3=1 > ( 17 ) 

where • • ■ denotes the measure part which we do not touch. The operator which 

is 

m— 1 



Kazakov introduced 10 



K = m 9r, 



m=l 



dl 



The action of K on the loop turns out to be 



■y oo £n oo Q 

Kw(£) = -j-H ~\ Y^m-*- — z n 



Use the Virasoro constraints eq. (Q) and undo the procedure to carry out the second 
derivatives of the couplings acting on the partition function. We find 

kwii) =« I" d£'tre £ ' M tre^ M » , (19) 
Jo 
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which becomes, after using the factorization property of singlet operators in the planar 
limit, 

Kw{£) = f w{£')w{£-£') . (20) 



This is Kazakov's loop equation ||10| . For more complete discussion as well as its 
extension to the supersymmetric loop equation and the determination of its critical 
points, see the recent preprint [[II] . 

Let us finally discuss how to take a double scaling limit of the Virasoro constraints 
at finite iV to obtain the Virasoro constraints indexed by half integers ( a single 
twisted boson) ||. First, the potential must be fine-tuned in order to reach the k- 
th multicritical point of Kazakov |10[ describing 2d gravity coupled to (2, 2k — 1) 
nonunitary minimal conformal matter. Let a be the lattice spacing in the level space 
of the orthogonal polynomials. As the leading contribution to the matrix elements of 
A is 2, we rescale as 

2-\ = a 2 / k \ sc , (21) 

with N — > oo and a — > 0, keeping 



9st 



1/ (a 2+l/k N) = finite . (22) 



In order to drive the system slightly away from the critical point, we add to the 
original potential a source which has a nontrivial continuum limit: 



£EM2-A 4 ) m/2 = £E^ _ (23) 



) m/2 = EE^ 

i I i i 

In the right hand side, we reexpanded the source in terms of the polynomial bases 
around the origin. This is added to the original couplings g^. We then undo the 
procedure to find a dressed ( renormalized ) source expressed in terms of the original 
bases {{(2 — \i) i+1 ^ 2 }}. Schematically, 



91 — > U + 91 — > 3i C) =U + --- ■ (24) 



We find 
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This is the Virasoro constraints of a twisted boson derived in [§] in the double 
scaling limit. It is worthwhile to emphasize that the limiting process can be postponed 
until the very end 0. 

B. Two Matrix Model 
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Let us now turn to the two matrix model. One starts with Mehta's formula ||18|| : 

dM 1 dM 2 e~ trvW( - Ml ' ) ~ trV< ' 2)(M2)+ctrMlM2 



Y[dXie *=i i A(Ai,...Ajv) (26) 



This time, we first discuss the approach based on orthogonal polynomials. Let us 
summarize some of the properties: 
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(27) 

(28) 
(29) 



Note that, unlike the case of the one matrix model, the number of terms appearing 
in the recursion relation (eq. fl29|) ) depends upon the degree of the potential. This is 
a crucial observation of |T3|, |T5| . ( See JT^, for earlier results.) 

It was Tada and Yamaguchi |14| ( and Tada pO| later) who first examined carefully 
the matrix elements of the Heisenberg algebra of the two matrix model: 



A, A 



c 



(30) 
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They were able to determine the potentials which can produce the critical points 
indexed by (p,p') = (4, 5), (3, 8), (3, 5), and (5,6), precluding explicitly the earlier 
expectation that two matrix model contains the critical points of the kind (3, *) 
alone. The same conclusion was reached by flT5f . 

Tada and Yamaguchi explicitly took the double scaling limit of the Heisenberg 
algebra: 



// 
A 

[Y,X] 



X 

Y = X. 
1 . 



q/p 



(31) 



which is distinct from the original more heuristic analysis of Douglas [pj. A prescrip- 
tion to obtain (p,p') critical points from an asymmetric potential of degree 2(p' — 1) 
has been given [BO]. A more general and detail analysis in the planar limit has been 



given recently in [pif] . 

Let us now discuss the approach based on the Dyson- Schwinger equation. It is a 
useful point of view to regard the second set of integrations in eq. ( p6|) as a Laplace 
transform £ of Ae 



_y(2). 
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(32) 



We again insert differential operators of degree m 



An (P) = £ f-^Y 

i (-ft \dfiij 



(33) 



Computing the action in two different ways, and using the formulas (eqs. (H),®) 
developed in the one matrix model, we find 
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w (n, m) 







(34) 



m + 1 



Res c=0 [C-Ql 



m+l 



31 + 



dji 



(-) 



n + l \ 



d 

32 + -7- 

d]2 



(35) 



(-) 



-9- 



Let me suggest how to take a double scaling limit of the constraint eq. (|34|),(|3~5|). To 
my knowledge, this procedure has not been carried out explicitly yet. Some progress 
has been made recently ||22|| . 

First of all, a set of couplings and gf^ must be found which produces the 
iPiP') critical point and we fine-tune the potential to this point. Insert a source term 
to the potential: 

- E {c-mY lP U- £ {c-\Y lP u . (36) 

(mod p) (mod p) 



Express the constraint eq. ( |3~3D in terms of the "renormalized" source in which the 
original couplings are absorbed. By taking the limit, we try to see the reduction of 
our constraints to the W p+ i constraints of |§ proven by |23 . 

Other results from loop equations include 
C. Discussion 

Here, we will discuss a few scattered items which deserve further investigations. In 
principle, one can apply the above method of deriving constraints to the c = 1 model. 
In practice, this is not such an easy task: we have to introduce a set of time dependent 
couplings to the potential of the one dimensional matrix model. This would lead to 
the functional differential constraints on the partition function even if we could write 
them down in a closed form. If we restrict our attention to time-independent ( static 



or zero momentum) quantities alone, it is possible to derive useful constraints [25 
The basic quantity of interest is the density of states 



P (m) = 4tr5 (m - h 



N 

1 00 / 1 \ 
= rrlm < n I x > < x I n > . (37) 

rather than the partition function. Or consider a once integrated quantity 

P (^) = — Li m tr In (// - h + id) . (38) 

The coordinate x of the one particle hamiltonian h originates from the eigenvalues of 
the matrix variable. Virasoro constraints are obtained from the reparametrization of 



eigenvalues but the procedure is not identical to that of the one matrix model [25 . 
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Our basic formulas (eqs. @, (|]), fl3l|), fl35|) ) can be regarded as a kind of bosoniza- 
tion of nonrelativistic fermions. A ramification to W geometry has been pursued 
recently [^6| . 

What we have been witnessing is the transmutation of the eigenvalue coordinates 
into the target space degrees of freedom. This phenomenon is also present in one 
dimensional matrix model [|27| , but conceptual understanding seems to be still lacking. 
Constraint equations in general provide nontrivial information on target space physics 
such as physics of black holes p8[ . Clearly much more work must be dedicated on 
these issues together with the issue of the definition of nonperturbative 2d gravity. 

I thank Joe Lykken for helpful discussion. 
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